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Abstract 

The intertwining operator constructed in |Szl[ ISz2| does not ap- 
pear in the right form. It is established there by using only the anticom- 
mutators J\ and J[. The correct operator involves all endomorphisms, 
J a , which are unified by the Z-Fourier transform. Although some of 
the correct elements of the previous constructions are kept, this idea 
is established by a new technique which yields the various isospectral- 
ity theorems stated in the papers on a much larger scale. The new 
results include new isospectrality examples living on sphere x ball- and 
sphere x sphere- type manifolds. Among them, there are such discrete 
isospectrality families where one of the members is homogeneous while 
the others are locally inhomogeneous (striking examples). Further- 
more, a large class of new isospectrality families are constructed by a 
deformations. 

1 Introduction. 

In papers [Szi~l lSz2j . the intertwining operator is constructed by the complex 
linear correspondence 

k* ■. tp(\x\,z)e Ql (x, Ji)...e Qp (A, Ji)e Qp+1 (A, j,)... e Qp+q (x, j x ) (i) 
- p(| n z)e Ql (x, j[) . . . q Qp (x, j[) . . . e Qp+1 (x, j[) . . . o Qp+q (x, 
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where Q±, . . . , Q p + q are arbitrary X- vectors, furthermore, @q(X, J\) = (Q + 
iJi(Q), X) and the corresponding Q'q(X, J{) are defined by the anticommu- 
tators Ji and J[ respectively. These anticommutators a-relate to each other, 
meaning, that the X-space is a direct sum, v = 

v (a) 0v (6) of k (a) 

resp. 

real-dimensional subspaces such that the components are invariant under the 
action of the anticommutators and they agree on the first component while 
Ji = — J[ holds on the second one. The complex linear property guaranties 
that it intertwines the operators d\D\» and diD[», which appear in the 
Laplacians involved to the angular momentum operators M = Yla=i 9 a D a » 
resp. M'. One of the ultimate goals is to intertwine the complete Laplacians. 

It was overlooked by the author that one can not allow arbitrary vectors, 
Qi, in the definition of k* because the operator becomes ill-defined. In fact, 
well-defined complex linear map can be introduced by choosing a system, 
{E^^E^}, where 1 < i < fcW /2; 1 < j < fc< 6 ) /2, of independent vectors 
in the corresponding component spaces v( a ) and which form a complex 
linear basis with respect to both complex structures J\ and J[. Then, the 
admissible Q's are the vectors laying in the real subspace spanned by all 
these E^'s. This natural complex linear map can be described in terms 
of the complex coordinates {z r } resp. {z' r } determined by the same basis 
{E r \l <r< + = k/2} = {E^} U {E^}, for the two complex 

structures respectively such that the image of a polynomial written up in 
terms of the coordinates {z r } is the polynomial of the same form but written 
up in terms of the other coordinates {z' r }. One can easily see that these well- 
defined maps depend on the real subspaces spanned by {E r }, meaning that 
maps defined for different real subspaces correspond different elements to 
the very same element in general. Thus, by allowing all possible Q's, the 
above k* is ill-defined indeed. 

Unfortunately, this problem can not be eliminated by replacing the ill- 
defined operator by one of these well-defined ones. In this case a much 
more serious difficulty appears, namely, such a well-defined operator does 
not intertwine d a D a » and d a D' a » satisfying a > 1, which are also parts of 
the corresponding Laplacians. Let it also be mentioned that the rest parts 
of the Laplacians as well as the boundary conditions are intertwined by it. 
In the papers, the proof of intertwining of the above parts of M resp. M' 
explores the false assumption claiming that the k* operates on functions 
defined for imaginary Q's in the same way as for the real ones. 

One of the reasons causing this blunder was that, instead of M and M', 
one was focusing just on D a * and D' a », i. e., tried to define intertwining 
operator using only the X-space. An other reason was that the intertwining 
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operator was defined just with endomorphisms J\ and J[. The corrected 
operator, involving all J a and focusing on M and M', is defined by choosing 
a basis {Ei, . . . , E k / 2 } and using Qj's laying in the real span of these basis 
vectors. More precisely, this correspondence is: 

k : f e [ ^ v ^(\X\,V)Q Ql ^X,V u )e Qp+1 ^ +q (X,V u )dV -» (2) 

J z 

[ ^ v (\X\,V)& Ql .jx,v u )@' Qp+1 p+q (X,V u )dV, 

Jz 

where &Q 1 ... P {X,V U ) := @ Ql (X,V u ) . . . Q Qp (X,V u ) and the corresponding 
0^ (X, V u ) are defined by the endomorphisms Jy u and belonging 
to the unit Z-vectors V u respectively. This operator associates functions 
defined by the Z-Fourier transform to each other. Appropriate intertwining, 
k~, can be established by using e~^ Z ' V ^ in the above formulas. Then also 
k r = (k + ~k~)/2 is going to be an intertwining operator. However, this 
paper proceeds only with the first version. 

The well-definedness of this operator follows from reasons such as only 
Q's laying in the real span of vectors Ei are used in its definition, further- 
more, the Fourier transform is an isometry on the corresponding 1? Hilbert 
spaces. The intertwining of the Laplacians is due to the Z-Fourier transform 
implemented into the formulas. The addition of this Z-Fourier transform to 
the original idea makes the mathematical situation much more complex, 
requiring a complete rethinking of the original construction. For instance, 
beyond proving the intertwining of the Laplacians, it is much more difficult 
to prove the intertwining of the boundary conditions for this operator. Ac- 
tually, the proof of the latter statement combines the Z-Fourier transform 
with an independent idea incorporated into the Independence Theorems. 
The ill-definedness of (pQ) was recognized by H. Fiirstenau whose observa- 
tion triggered the author's thorough rethinking of his complete construction. 
The much deeper problems hidden under the cover of ill-definedness came 
to the light during this revision process. The reborn operator presented here 
saves all the previous results and provides also new interesting isospectrality 
examples. The main goal is to establish 

Theorem 1.1. The ball- resp. sphere-type manifolds, which have the same 
radius function ip(\X\, \Z\) and are defined on H-type groups H^ b) having 
the same parameters a + b and I, are isospectral. This statement extends to 
a large class of general 2-step nilpotent Lie groups where an isospectrality 
family is defined by a -deformations. 
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On ball-type manifolds this statement includes the well-definedness of ^ 
and the intertwining of the Laplacians as well as the boundary conditions. 
Since the Dirichlet condition is intertwined, by restrictions, the operator 
induces bisections between the function spaces defined on the boundaries. 
By observing that it is enough to use functions satisfying the Z-Neumann 
condition on the ambient manifolds, one can prove the intertwining property 
also on the boundary manifolds. 

All these statements extend onto the solvable extensions of 2-step nilpo- 
tent Lie groups. Furthermore, new examples, not discussed in the original 
papers, are also constructed. They live on spherexball- and spherex sphere- 
type manifolds. Among them two are particularly interesting. Namely, the 
isospectrality family of spherex sphere-type manifolds, constructed both on 
H^ a ' b ^ and SH^ a ' b \ the metric is homogeneous for the manifold belonging 
to the pair (a + 6,0) or (0, a + b), while the metrics satisfying ab / are 
locally inhomogeneous. Also the dimensions of groups of isometrics acting 
on the members are different. These are new contributions to the old list of 
striking examples constructed on the spherex torus-type manifolds of ifg . 
resp. geodesic spheres of SH^ a ' b ^ . 

These theorems are established, first, on H-type groups and their solvable 
extensions. Then, they are extended to those 2-step nilpotent groups and 
their solvable extensions which are defined by endomorphism spaces obtained 
by perturbing the endomorphism space of a given H-type group . 

The perturbation process mentioned above is as follows. The endomor- 
phisms, Jz, on H-type groups are defined by endomorphisms, jz, acting on 
the irreducible components, by the formula Jz = {jz, ■ ■ ■ ,jzi —jz, ■ ■ ■ , —jz)- 
The perturbation primarily concerns the endomorphisms jz, i- e., close to 
the Cliffordian one, a new linear space of endomorphisms is chosen with ele- 
ments denoted by jz- Then, the groups, H\ a,b ^ and SH\ a ' b \ resulted by such 
a perturbation arise from endomorphisms Jz defined by the same formula 
in terms of jz- 

Such a perturbation results a new family, defined by the same (a + b) 
and I, whose members are obviously a-equivalent. By a theorem proved in 
the last section, the k always intertwines the Laplacians, however, the same 
statement for the boundary conditions is not guarantied. This problem is 
solved by the independence theorems which state the independence of certain 
subspaces formed by functions. This independence is established, originally, 
for H-type groups defined by Cliffordian endomorphism spaces and remains 
true for groups defined by endomorphism spaces which are close enough to 
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the Cliffordian ones. The latter spaces constitute an open set whose mem- 
bers are called the small perturbations of a given Cliffordian endomorphism 
space. The complete isospectrality theorems are established only on those 
a-equivalent groups whose endomorphism spaces are produced by small per- 
turbations of the Cliffordian ones. It is also important to understand that 
the independence theorems alone do not validate the intertwining property 
for the boundary conditions. They can guarantee it just for intertwining 
operators defined by cr-deformations. 

The new features in this mathematical process include, first of all, the 
integral formula by which the intertwining operator is defined. This formula 
deeply roots in quantum theory [Sz4[ ISz5] . furthermore, it can be used also 
for explicit computations of the eigenfunctions and spectra. This rooting 
in physics is exhibited by the surprising fact that the Laplacian on the in- 
vestigated manifolds can be identified with the Landau-Zeeman operator 
attached to electron-positron systems where these particles are orbiting in 
constant magnetic fields. Above, the endomorphisms — jz resp. jz corre- 
spond to electrons resp. positrons, and cr-deformations is interpreted such 
that some of the electrons are exchanged for positrons. Yet, the spectra 
on all submanifolds investigated in this paper are not changing during this 
exchange-process. The local geometry, however, is dramatically changing. 
The manifolds in the striking examples, for instance, are homogeneous for 
systems having particles of the same type, while, they are locally inhomo- 
geneous for mixed particles. 

This is a physical interpretation of the above isospectrality theorems. 
The perturbation can be interpreted such that, instead of a system of iden- 
tically charged particles, one considers ones which are charged distinctly. 
Then, the isospectrality theorems are established also for systems produced 
by small perturbation of the charge. The fact of non-changing spectra during 
electron-positron-exchanges is well known in physics. However, the state- 
ment of this form concerns the spectra considered on a non-compact mani- 
fold. Our statement claims much more than just this. Namely, the spectra 
remain the same also on a large class of compact submanifolds. An other 
distinguishing feature is that no attached local geometries are considered in 
quantum theory. 

The intertwining of the Laplacians can be established by using only the 
integral formula. For proving the intertwining property for the boundary 
conditions, a new idea, appearing in the Independence Theorems, is involved. 
These theorems are also important new features in this field. 

The methods developed in this paper apply only for cr-deformations. 
A characteristic feature of these discrete deformations is that they do not 
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change the Ricci curvature. This fact is strongly used in proving the in- 
tertwining regarding the boundary conditions. This experience strongly in- 
dicates that the submanifolds considered in this paper are not isospectral 
on the Gordon- Wilson [GW| examples where continuous isospectral defor- 
mations with changing Ricci tensor are established on groups defined by 
2-dimensional Z-spaces. Other arguments supporting this statement are ex- 
plained in the end of the last section. 

The much more general Theorem 11.11 replaces the Isospectrality The- 
orems of the articles. Only the construction of the intertwining operator 
(cf. pages 461-465 in [Szl] and 371-375 in [Sz2j ) is effected by this problem. 
The major non-effected part includes all the Non-Isometry resp. Rigidity 
Theorems and the preparatory part of Sections 4. resp. 3.. This problem 
with solution was announced at the CUNY Geometric Analysis Conference, 
in 2006 [5z3] . 

2 Technicalities. 

The constructions are performed on 2-step nilpotent metric Lie groups and 
their solvable extensions. The nilpotent groups are defined on correspond- 
ing orthogonal direct sums, v © z, of Euclidean spaces where the compo- 
nents, v = R k and z = R' are called X- and Z-space respectively. The 
Lie algebra is completely determined by the linear space, J z , of skew en- 
domorphisms whose actions on the X-space are defined by the relation 
([X, Y],Z) = (Jz(X),Y) , where X, Y £ v and Jz is the endomorphism 
associated with Z £ z. The Riemannian metric, g, is the left invariant ex- 
tension of the natural Euclidean metric on the Lie algebra. The exponential 
map identifies the Lie algebra with the group itself, thus also the group 
can be considered to be defined on the same (X, Z)-space. Each group, 
(N,g), extends into a solvable group (SN,g s ), where a point is represented 
by (t,X,Z). 

Particular 2-step nilpotent Lie groups are the so called Heisenberg-type 
groups, defined by endomorphisms satisfying the Clifford condition J| = 
— \Z\ 2 id. These metric groups are attached to Clifford modules, thus the 
classification of these modules provides classification also for the H-type 
groups. In this case the X-space decomposes into the product v = (R r 'W) a + 6 : 
R r W a x R r W fc and endomorphisms Jz are defined by endomorphisms jz 
acting on the smaller space R r (^ such that they act on R r W a resp. R r W fe 
according to the Cartesian product jz x • • • x jz resp. —jz x • • • x —jz- 
The H-type groups are denoted by the symbol H\ a,b \ which indicates the 
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above decomposition. The solvable extensions of H-type groups are denoted 
by^ (a ' b) . ' 

The Laplacian on a H-type group is of the form 

1 r 
A = A x + (l + -\X\ 2 )A z + Y,d a D a ; (3) 

a=l 

where D a » denotes directional derivative along the vector field X — ► J a (X) = 
Jz a (X) defined for each element, Z a , of an orthonormal basis of the Z-space. 
The integral curves of this field are the Hopf circles defined for the complex 
structure J a . In the isospectrality constructions performed in this paper one 
should deal with this compound operator. Earlier, the constructions were 
performed on center periodic H-type groups, T\H, defined by factorizing the 
center of the group with a Z-lattice r = {Z 7 }. In this case the L 2 function 
space is the direct sum of function spaces W 7 spanned by functions of the 
form ^ 7 (X, Z) = ip(X)e 2nl ( z ~" z ' ) . Each W 1 is invariant under the action of 
the Laplacian, i. e., A^ 7 (X, Z) = a y ip(X)e 27T ^ z -" z ' } , where operator D 7 , 
acting on L 2 (v), is of the form 

□ 7 = A x + 2ttLD 7 • -4vr 2 |Z 7 | 2 (l + j\X\ 2 ). (4) 

Note that the first operator involves all endomorphisms J a while the second 
one involves, regarding each invariant subspace W 7 , only J 7 . 
Remark. There is pointed out in |Sz3|, ISz41 ISz5] that operator @ is noth- 
ing but the Zeeman-Hamilton operator of a free charged particle (the 2D 
version is called Landau Hamiltonian) , which was used for explaining the 
Zeeman effect. Term involving L> 7 » is the so called angular momentum 
operator, which represents a preliminary version of the spin concept. The 
non-periodic metric group (N, g) strongly relates to Dirac's relativistic multi- 
time model, which, in order to furnish relativistic features on the quantum 
level, endowed the particles with individual self-times. In the H-type model 
the multi-time is represented by the multi-dimensional center of the group. 
Regarding this relativistic interpretation, which is not the same as the clas- 
sical relativism, Laplacian ([3]) on the total space (space-time) corresponds 
to the Klein-Gordon operator. Note, however, that this multi-time operator 
is an elliptic one. This fact points to the distinctive features of the multi- 
time and classical relativism. Operator, M, involving all angular momentum 
operators is called compound angular momentum operator. 
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3 Isospectrality constructions. 



The isospectrality constructions are performed on H-type groups, H\ ' , 

and on their solvable extensions, SH^ a,b \ first. It is only the last section 
where they are extended to u-equivalent groups whose endomorphism spaces 
are constructed by perturbing the Clifford endomorphism spaces. The main 
goal is to describe these constructions on non-periodic groups, however, in 
order to see both the similarities and differences, they are briefly reviewed 
here also in the center periodic cases. For fixed a + b and I, these groups 
are defined on the same (X, Z)- resp. (i, X, Z)-space. There is established 
in many different ways that the metrics, gf' b \ in a family have completely 
different local geometries (cf., for instance, the striking examples), yet they 
are isospectral on a wide range of submanifolds. 

3.1 Constructing the ballxtorus- and sphere x torus-examples. 

These examples are constructed for a family, T\H^ a,b \ of Z-periodic man- 
ifolds. The submanifolds considered are torus bundles over a ball (resp. 
sphere) around the origin of the X-space. An intertwining operator can 
be constructed such that, for each invariant space W-y constructed above, 
just an orthogonal transformation conjugating J 7 a '^ to J 7 a ' b ^ on the X- 
space should be considered. The intertwining operator on W 7 is defined 
by the map induced on functions ip(X), defined in the Fourier- Weierstrass 
decomposition, by this point transformation. This transformation clearly in- 
tertwines d 7 with □' such that it keeps also the boundary conditions. (The 
boundary conditions can be described in terms of radial functions. The 
intertwining of boundary conditions is then due to the invariance of these 
functions under the action of the operator.) It induces an appropriate in- 
tertwining operator also on the boundary manifolds. The striking examples 
appear on the quaternionic families H^ a ' b \ in which case the sphere x torus- 
type boundary manifolds in T\H^ a+b ' ^ are homogeneous while the others 
in the family are locally inhomogeneous. Note that the simplicity of this 
case is due to the fact that the intertwining operator is constructed, on each 
invariant space W 7 separately, by a single endomorphism, J 7 . 

3.2 The ball- and sphere-type domains. 

These examples were originally constructed in [Szll Sz2j. The ball-type 
domains are, by definition, diffeomorphic to Euclidean balls such that the 
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sphere-type boundary manifolds are level sets described by equations of the 
form \Z\) = resp. (/?(|X|, \Z\, t) = 0. These domains are invariant 

under the action of the orthogonal group 0(R*) x O(R'), thus, they may 
be called domains of (X, Z)-revolutions. They can be visualized such that 
there is an X-ball centered at the origin of the X-space considered over the 
points of which there are Z-balls of radius around the origin of 

the Z-space considered. Then, the boundary is a level set defined by the 
equation <^(|X|,|Z|) = \Z\ — Rz(\X\) = 0. By this reason, function (p is 
called radius function. 

Note that radius is constant along a sphere defined by a constant 

radius Rx = \X\ in the X-space. Furthermore, the ball bundle defined by 
the Z-balls over this X-sphere is trivial. These are the so called spherex ball- 
type manifolds whose boundaries are spherex sphere-type manifolds. The 
new examples, not discussed in the earlier papers, are constructed on these 
domains and surfaces. 

An other visualization can be started out with a Z-ball in the Z-space 
over the points of which there are X-balls of radius Rx(\Z\) considered. 
Then the boundary is defined by |X| — _Rx-(|Z|) = 0- However, this paper 
proceeds with the first description. 

In the solvable case one should consider (Z, i)-balls and (Z, t)-spheres 
around the origin (0^,1) defined for the hyperbolic (Z, t)-space. The base 
manifold is the same X-sphere as before. Note that a Z-ball Br z (0z) (resp. 
Z-sphere Sr z (Qz)) uniquely extends into a geodesic ball (resp. sphere) of 
the hyperbolic (Z, t)-space. A spherex ball- type domain can be described as 
a hypersurface in a ball-type domain such that the Z-balls (resp. (Z, t)-balls) 
of the ball type domain are considered only over the points of a sphere Sr x 
laying in the X-space. Similarly, the spherex sphere- type manifolds can be 
regarded as hypersurfaces in the sphere- type manifolds. 

The isospectrality will be investigated, first, for the discrete families, 
ff( a ' b ) ^ defined by the same a + b and I. The Laplacian is described then by 
(|3|). Comparing with the Zeeman operator (JH), this operator involves all the 
endomorphisms, making the constructions much more difficult. The Lapla- 
cians of the members in a family differ from each other just by the last term, 
M, which is called compound angular momentum operator. The spectral 
investigation both of M and A is completely missing in the literature. Note 
that this most intriguing operator, M, commutes with both operators in the 
rest part of ([3]). 
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3.3 Eigenfunctions motivating the intertwining operators. 

The eigenfunctions constructed next are not directly used in the isospectral- 
ity constructions and the rest part of this paper is understandable without 
knowing about their actual explicit form described in the second half of 
this section. However, there are important concepts introduced in the first 
part which are heavily used later on. The ultimate reason for describing 
these functions here is that they very clearly suggest the explicit form of the 
sought intertwining operators. 

Since M commutes with the rest part, O, of A, the eigenfunctions of A 
can be sought as common eigenfunctions for both operators M and O. In the 
very first step we look for the eigenfunctions of a single angular momentum 
operator Dy», defined by a Z-vector V. For a fixed X-vector Q and unit 
Z-vector V u = j^V, consider the X-function Qq(X,V u ) = (Q + iJ\/ u (Q), X) 
and its conjugate Qq(X, V u ). For vector V = \V\V U , these functions are 
eigenfunctions of Dy with eigenvalue — |V|i resp. |V|i. The higher order 
eigenfunctions are of the form ©q©q with eigenvalue (q — p)\V\i. 

In order to find eigenfunctions of the compound operator M, consider a 
sphere Sr z of radius Rz around the origin in the Z-space. For an appropriate 
function 0(|X|, V), depending on \X\ and V 6 Sr z , define 

FQ pqRz (cf>)(X,Z)= <f e^ z ' v ^{\X\,V)%{X,V u )Q q Q {X,V u )dV. (5) 
Js Rz 

By V u = V/\V\, the V u is considered as a function depending on V. Due to 
the relation M <f = <f iDy, this function is an eigenfunction of M with the 
real eigenvalue (p—q)Rz- These functions are eigenfunctions also of Az with 
eigenvalue R z . Also note that these eigenvalues do not change by varying 
Q, or, if the simple functions Qq(X, V u ) resp. Qq(X, V u ) are exchanged for 
their pluralistic versions @Q 1 ... P (X,V U ) := &q 1 (X, V u ) . . . Qq p (X, V u ) resp. 
®Q p+ i... p+q (X,V u ). 

Functions © defined by simple resp. pluralistic functions are said to 
be one-pole resp. multiple-pole functions with poles Q resp. {Qi}- The 
function space generated for fixed 1-pole (resp. multi-poles) by all possible 
<j) is not invariant with respect to the action of Ax, thus the eigenfunctions 
of the complete operator A do not appear in this form. In order to find the 
common eigenfunctions, the homogeneous but non-harmonic 1-pole polyno- 
mials QqGq (resp. the multiple-pole polynomials) of the X-variable should 

be exchanged for the 1-pole harmonic polynomials I1x\@q@q), (resp. to 
the corresponding multiple-pole harmonic polynomials) defined by the pro- 
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jection, Ilj^ , onto the space of n = (p + g)-order homogeneous harmonic 
polynomials of the X-variable. These projections are explicitly described in 
the form 

=A° X + B[ n) \X\ 2 A x + B^ n) \X\ 4 A x + ... (6) 

in [Sz2| (cf. formula (3.14) there), where A^ = id and the constants B^ 
are determined by a recursion formula. This formula easily implies that also 

KF QpqR M{X,Z) = I S zy U{\X\,V)n x (%{X,V u )Q q Q {X,V u ))dV 
Js Rz 

(7) 

are eigenfunctions of M and Az with eigenvalues belonging to ([5]). The 
same statement is true regarding the multiple-pole-cases. 

The action of the complete Laplacian (|3]) is a combination of X-radial 
differentiation, d\ x \, and multiplications with functions depending just on 
|X|. I. e., the action is completely reduced to X-radial functions. Also this 
reduced form of the Laplacian is not changing by varying Q, or, switching to 
multiple-pole functions. The eigenfunctions of A can be found by seeking the 
eigenfunctions of the reduced operator among the X-radial functions. The 
explicit computations are carried out in |Sz4|, ISz5] . Since these details are 
not used in this paper, we just indicate that the eigenfunctions appear in the 
form f s e^ z ^ (p(V)F^' q) (X,V u ))dV, (resp. in a corresponding multiple- 
pole version of this function), where Fq' q \x, V u )) is an eigenfunction of 
operator D 7 satisfying |F 7 | = |V| = Rz- The latter ones are explicitly 
described in [Sz41ISz5] in terms of homogeneous harmonic polynomials which 
are multiplied with radial functions. The eigenfunctions are determined 
below also by a different method, using the so called ltd polynomials. 

Note that this construction is carried out for a fixed 1-pole Q (resp. a 
fixed multiple-pole, {Qi, . . . , Q p , Q p +i, • • • , Q p + q })- An other type of con- 
structions is as follows. For any unit vector V u of the Z-space, consider a 
complex orthonormal basis {Qv n i, ■ ■ ■ > Qvuk/2} on the complex X-space de- 
fined by the complex structure Jy u such that the vectors in front lay in v( a ) 
and all the others are in \~( b \ Such a basis defines the complex coordinate 
system {z Vul = &q 1Vu , . . . , z Vuk/2 = &Q (k/2)v J on the X-space. This ba- 
sis field must be smooth on an everywhere dense open subset of the unit 
Z-sphere such that it is the complement of a set of measure. For given 
values p\ , qi , . . . , , Qk/2 1 consider the polynomial nf=i z v i^V i • Then the 
functions 

k /2 

I e^UiVMlXl) fl^ ui 4 ui dV (8) 
Js R Z i=l 



11 



are eigenfunctions of the compound angular momentum operator M. In 
order to have an eigenfunction for the complete Laplacian, one can use 
the above described method of projecting the polynomial into the space of 
homogeneous harmonic polynomials which have order {p\ + q\ + • • • +Pk/2 + 
■ In |Sz4t ISz5] . the eigenfunctions of D 7 are determined also by an other 

method, seeking them in the form hy iq ^ e~ Rz ^ ZVu,Zv ^, where the hy l<1 ^ is 
an (pi + qi + • • • + p k / 2 + (/fc/2)-order polynomial. Then this function is an 
eigenfunction of D~ satisfying \Z~\= Rz if and only if the latter function is 
an ltd polynomial regarding the complex structure Jy u . The final form of 
the eigenfunction is 

<f S zy) <t>{V)h^ qi) e- RzZv ** v *dV. (9) 
Js Rz 

Since Ito's polynomials are non-homogeneous, this is a different representa- 
tion of the eigenfunctions. These explicit descriptions of the eigenfunctions 
will not be used in the rest part of the paper. 



4 Constructing the intertwining operators. 
4.1 Constructions on ball-type domains. 

The constructions described in this sections are carried out for Heisenberg 
type groups, H^ a ' b ' and H^ a ' b \ which are in the same isospectrality family, 
i.e., a + b = a' + b' holds. From each eigenfunction-construction described 
above one can derive the corresponding intertwining operator intertwining 
the corresponding eigenfunctions provided by the construction. Note that 
the functions appearing there are not of class I? regarding the Z-variable 
thus integral formulas (jS} , Q , ([H} can not be directly used for defining the 
operator. This is why function 0(|X|,V), depending on \X\ and V £ Sr z 
originally, is exchanged for one which an L 2 function of the V-variable, for 
any fixed \X\, and the integral is taken over the whole Z-space R'. In other 
words, the Z-Fourier transform on L^-setting is considered. 

Also the order for introducing the various versions of the intertwining 
operators is an important issue. The first version is defined for a fixed basis, 
Qf = {Qii ■ ■ ■ 3 Qk/2} = {Qj? 3 Qi^}) wn i cn does not depend on V u , where 
the first /2 number of vectors are in v( a ) and the following /2 number 
of vectors are in \~( b \ If these vectors are chosen such that they form an 
orthonormal basis regarding a fixed Jy 0u , then they form a complex (in 
general, non-orthonormal) basis for V u 's which form an everywhere dense 
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open subset on the unit sphere of the Z-space. This operator is defined by 
means of the polynomials written up in terms of the coordinate functions 



z Vul (X) = @ Ql (X,V u ),...,z Vuk/2 (X) = Q Qk/2 (X,V u ) 



(10) 



resp. {z' Vul , ■ ■ ■ , z' v k / 2 }- The denotation indicates that, although the basis 
is fixed, these coordinate functions depend on V u . 

For constructing the eigenfunctions in the previous section, one is using 
a fixed Z-sphere and functions, (j)(V), of Dirac-type concentrated on this 
sphere. Then the eigenfunctions are represented by Z-Fourier transforms of 
these Dirac-type functions. Whereas, in defining the intertwining operator 
kq f determined for a constant basis Q_p, one considers Z-Fourier transforms 
of appropriate L^-functions. More precisely, for an L|-function, 4>(V), de- 
fined on the Z-space and arbitrary set {pt,qi}, where i = 1, . . . ,k/2, of 
natural numbers the intertwining operator is defined by 



I. e., the kq f corresponds to a function, which is defined by the Z-Fourier 
transform formula in terms of (f), zf l , ~zf , the function defined by the same 
expression but which is written up in terms of cj), z'f 1 , and zf* . In these for- 
mulas, the V u = V/\V\ is considered as a function depending on V, further- 
more, the dependence of the complex coordinate functions on the X-variable 
is described in (llOp . Note that in this first version of the operator function 
4> depends just on V and not on \X\. 

The domain and range of this operator is discussed in the next section. 
In this section one is focusing on the well-definedness and the intertwining 
of the Laplacians. Concerning these questions, we have. 

Theorem 4.1. The above kq f is a well defined one-to-one operator. 

Proof. This theorem is well established by proving that the image of a func- 
tion which is in the domain of kq f and vanishes almost everywhere is a 
function vanishing almost everywhere. For proving this statement, suppose 
that function 



where the terms of the sum (series) are defined regarding the independent 
polynomials Yii=i z v u i^v u ii van i snes almost everywhere. Since the Z-Fourier 




k/2 
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transform is a one-to-one map on the corresponding L 2 -Hilbert space, for 

any fixed X, function ip(X, V) = Z{ Pi , qi} hnm}^) UiH %JX)1*JX) 
(whose Fourier transform is considered) must vanish almost everywhere. 

By the independence of the polynomials Yli=i z v i^v i> wria t is satisfied for 
almost all V u , a general function <p defined by this formula is non-zero if 
there is a non-zero L 2 -function 4>{pi,qi} among the component functions. 
Therefore, due to the assumption of this theorem, all these </>'s must vanish 
almost everywhere. Then, also the image tip' must vanish almost everywhere. 
This proves the statement completely. 

Actually, one has proved the following stronger statement: Function 
tip' is zero almost everywhere if and only if its preimage tip is zero almost 
everywhere. Thus also the one-to-one property is established completely. □ 

Much more handy alternative definitions of the very same acq f are es- 
tablished in the following theorem. 

Theorem 4.2. (A) In the definition of kq f , function (ft may depend also on 
\X\, or, even more, on \X^\ and \X^\, where X = X^+X^ corresponds 
to the decomposition v = v( a ) © . I. e., an equivalent version is: 

fc/2 



KQp :f QF{Mi} W(I,Z)= f J< z M4(\X\,V)l[%jX)?%JX)dV (12) 

1 i= i 

-^MwtfttX*.*) = / ^ y ^{\X\,V)^4\lX)il^X)dV. 

i=l 

In the more general version, the </>(\X\,V) is replaced by 0(\X^\,\X^\,V). 

(B) The kq f intertwines both the Euclidean Laplacian Ax and the pro- 
jections Tlx , n^" a \ Ilj^ with themselfs respectively. These projections are 
described in furthermore, n = p + q = J^Pi + Y^Qi> n a = Pa + Qa = 
J2i=i(Pi + 9i), n b =p b + q b = J2i=a+i(.Pi + 9i)- The operator can be written 
in the following alternative form: 

k /2 

K QF :-HF QF{mi} mX,Z)= f e i ^ v U(\XlV)U%\l[z^JdV (13) 

JRl 1=1 



In the more general version, the 4>(\X\, V) resp. n^(fJi=i • • • ) are replaced 

Ia+b \ 
1=0+1 ■ ■ ■ )■ 



by ^(\X(%\X(%V) resp. U { ^\UU ■ ■ ■ )^x\Ui +b 
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(C) These versions for defining kq f allow to introduce its domain in a 
more precise way. In order to work on I? function spaces, one should con- 
sider functions of the form <fi(\X\,V) = e~' x ' where (p, for any 
fixed \X\, is of class I? with respect to the V -variable, and it is a polynomial 
with respect to the \X\ -variable. By plugging them into the Fourier trans- 
form formula, they generate a pre-Hilbert space whose closure, regarding the 
l?-Eilbert norm, is a Hilbert space. (In the next section, this domain is 
identified with the standard I? -Hilbert space defined on M. h A larger 
domain can be generated by functions (f>(\X\,V), where, keeping the above 
assumption regarding the V -variable, function c/)(\X\,V) depending on vari- 
able \X\ is of class I? for almost all fixed V on any interval < \X\ < R. 

Proof. (A) This proof explores that the system Qp of independent vectors 
decomposes into two subsystems, Qf = {Q*- a \ Q^- 1 }, consisting vectors from 
v*- a ) resp. v^ b \ They form a complex basis, both for Jy u and Jy^, for an 
everywhere dense open set of unit vectors V u . This basis can not be orthonor- 
mal for all V u , even though it is orthonormal for some V u . For each V u , let 
Ry u = {R^,Ry^} be a complex orthonormal basis, regarding the complex 
structure Jy u , defining the complex coordinate system {zr v j}. Then, there 
exist complex matrix, (cjj(V u )), such that z-R Vu i = J2j c ij (^u) ZQv u j hold. 
The equation expressing the orthonormality is ^^ /mn Ci rn Cj n ( i Q rn ,Q n ) — 
where matrix with entries {Q m , Q n ) is real. The basis field can be chosen 
such that it is continuous on an everywhere dense open subset of the unit 
vectors V u . 

In terms of J v , the above matrix transformation defines an other com- 
plex coordinate system, z' R/ i = ^ - (Hj(V u )zQ F y;, with the corresponding 

complex basis Ry . Although it is not the same as Ry u , but, due to the 

relations Jy^ = Jy^ and j'y^ = —Jy\ this new basis also yields the above 
orthonormality equation. Therefore, it is orthonormal regarding the com- 
plex structure J' Vu and, according to the computations: 

\X\ 2 = J2zR Vu iZR Vu i = ^(^(HjiV^ZQpV^^CiriV^ZQ^r) = (14) 

i i j r 

= C ij( V u)Cir{V u ))zQ F V u jZQ F V u r) 
j,r i 

= ^2 Z Rv u i^Rv u i = ^(^ C ij( V u) Z ir(V u ))ZQ FVuj ^Q FVur ) 
i j,r i 
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the X-radial function \X\ 2 appears in the same polynomial form regarding 
both coordinate systems {zQ F v u j} and { z Q F v u j}- ^he same statement is 
true for any power, or, by using power series, for any function </>(|X|) of the 
basic radial function. This observation proves (A) completely. 

(B) Regarding the Laplacian the same computation yield: 

i i j r 

= EE c «j(K)c> r (K))9 Z[if ,„A qpV „ r ) 

j,r i 

=zZ d ^ Vu M Vui = E(E c ^(^)^(^))^ f v 11 Aq^ 11 .)' 

i j,r i 

i. e., also this Laplacian appears regarding both coordinate system in the 
same form. This proves the invariance of Ax under the action of kq f . The 
explicit formula §U§ regarding ITj^ along with (A) and the above statement 
concerning Ax prove (B) completely. 

(C) This statement is self-contained. □ 

Yet an other alternative definition of kq f can be introduced by using 
1-pole functions 0q0q, where Q is in the real span of the vector-system 
Qf, i- e., Q G Span^Qp). The version using multi-pole functions, defined 
by vectors laying in Span^Qp), adds nothing new to the above polynomial- 
version, therefore, this case is omitted here. Also note that the multi-pole 
functions span the same function space spanned by the 1-pole functions, 
thus, complete analysis can be performed by using only the simple ones. 

Theorem 4.3. The kq f is well defined by each of the following correspon- 
dences: 

KQ F : Fq p M(X, Z) = f Rl e^y)^V)% Vu {X)Q q QVu {X)dV (16) 
- F^mX, Z) = e^y)^V)Q% Vu (X)Q% Vu (X)dV, 

*Q F ■■ F Qpq {MX,Z) = f Rl e^ v )ci>(\X\,V)% Vu (X)e q QVu (X)dV (17) 

- r Qm mx,z) = f Rl e ^ v )H\xiv)e^ Vu (x)e' q Vu (x)dv, 

*Q F ■ Wa«(0(*, Z) = f Rl e^y)miV)Ii% +q \% v Q q QVu )dV (18) 
- UT' Qm mX Z) = f Rl e^V)m\,V)^\®% v ®'« Vu )dV, 
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where Q 6 Span^(Qp) is an arbitrary vector and Qqv u (X) := Qq(X,V u ). 

There is a more general version also in this case which corresponds to the 
exchange of Q for a pair, (Q^ a \Q^), followed by the exchange of functions 
<j>(\X\,V), % Vu %V u > n^ +9) (..) for the following ones 

m { % \x®\,v), %W^W„%V u ^V a , nfer J (..)nfer } (-) 

respectively. 

Remark. Although, it is not defined by an appropriate basis, the well- 
definedness of the operator is not jeopardized in this theorem. It can be 
defined, however, by constructing a basis as follows. 

For fixed Q and natural numbers p and q, let *&Q pq (resp. Hq P9 ) be the 
I? function space spanned by functions of the form J r Q pq ((j))(X, Z) (resp. 
T-CJ-Qpq(<p)(X, Z)), where 0(|X|,V) can be an arbitrary L^-function. For 
fixed p and q, all these spaces sum up to the total spaces <& pq = J2q ^Qpq 
(resp. S pq = SQpg). There exist finite many Qi such that the total 
space is the direct sum of the independent subspaces *&Q m (resp. 3q iP(J ). 
For the two type of total spaces these numbers are different. Due to the 
non-degeneracy of kq f , the total spaces <f>' pq (resp. E' pq ) are the direct sums 
of the independent subspaces &Q iPq (resp. S'q ). The kq f can be defined 
just by its actions k Qf : ® QiPq & Q . pq (resp. k Qf : S QiPq 5' Q . pg ) on 
these subspaces. This construction method, which will not be used in this 
paper, can be applied for explicit spectral computations. □ 

Above, six versions of the very same intertwining operator defined by 
a constant complex basis, Q^, were introduced. A changing orthonormal 
complex basis field, Q(V U ) = {Q {a) (K), Q {,,) (K)}, which is supposed to be 
continuous on an everywhere dense open subset of the unit vectors V u , also 
defines an intertwining operator. Unlike the constant field, which can be 
orthonormal only for V u 's of zero measure, the changing field is supposed 
to be orthonormal almost everywhere. Then, by Jy^ = Jy^ and Jy = 

—Jy \ also the basis Q,'(V U ) is orthonormal regarding the complex structure 
J'y^. The starting version of the intertwining operator k-q(v u ) defined by a 
changing complex orthonormal basis is introduced by formula pip , where 
the complex coordinates, zy u i, are defined by Q(V U ). Now we have 

Theorem 4.4. Operator 

K Q(V U ) ■ FQ(V u ){ Pi , qi }(<t>)(X, Z) -> ^Q(y u) { Pi , gi }(<£)(^, Z), (19) 

defined for all sets, {pi,qi}, of natural numbers and L 2 -functions <j>(V), is 
a well-defined one-to-one map leaving the X-radial functions, the Laplacian 
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Ax and the projections U x invariant. Thus, this operator can be defined 
in the alternative ways 

K Q(Vu) ■ ^Q(K){ P< , gi }(^)(^ Z) -► ^Q( Vu ){ Pi , qi }^)( X ^ Z )i ( 2 °) 
K Q(v u ) ■ KFQ(v u ){ Pi , qi }((l>)(X,Z) ^HJ r, Q{Vu){nm} {4>){X,Z), (21) 
where the <f>(\X\,V) may depend also on \X\, or, \X^\ and \X^\. 

The proof is the same as for the constant basis case. Since both basis', 
Q(V^) and Q'(V U ), are orthonormal, the proof of the invariance of the radial 
functions, Euclidean Laplacian Ax, and projections is even simpler as 
in the previous case. Let it also be mentioned that versions (|16|) - f)18j) can 
not be introduced for the changing basis case because the Q's must be in 
the intersection of all real subspaces span^Q(V u ). 

Now we are ready to prove the first main theorem in this paper. 

Theorem 4.5. Both in the constant and the changing basis cases, the kq 
intertwines the complete Laplacians A and A'. 

Proof. There is proved above that both the Laplacian Ax and operators de- 
fined by multiplication with radial functions are intertwined by the kq . The 
other parts of the Laplacians are also intertwined because of the following 
identities. 

k q : M^ Q(K){p . )(? . } (</>) = ^Q{ Pi)9i }((g - p)\V\<f>) (22) 
k q : A x JP Qfei3i} (0) = F Q{Pi , qi} (-\V\ 2 (l>) (23) 

where p = ^2pi and q = Ylli- These formulas remain true if T is exchanged 
for HT. □ 

4.2 Constructions on sphere x ball-type domains. 

For introducing the intertwining operators on sphere x ball- type domains, 
one can start with version (|13|> . where it is defined for a constant basis 
in terms of homogeneous harmonic polynomials of the X-variable defined 
on the ambient space. It is pointed out there, that the operator is well 
defined without using a basis, but now, for each n, consider functions of 
the form n£ \Ui=i z v u SvJ = \ X \ 2n a v u { Pl , qi } such that, for a fixed V u , 
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functions <?v u {Pi,qi} f° rm a basis among the corresponding spherical har- 
monics defined on the unit sphere of the X-space. Note that the dimen- 
sion of n t?l -order harmonic polynomials is less than the dimension of n th - 
order homogeneous polynomials, thus, not all polynomials from the lat- 
ter set are subjugated to the projection. Anyhow, such choices for such 
basis' exist. All those V^'s satisfying this property form an everywhere 
dense open subset of the unit Z-vectors. The functions whose Fourier trans- 
forms are considered in an arbitrary version of the intertwining operator 

have unique expansions, Efe,^} ( t ) {p l ,q 1 }(\ x \^ v ) cr v u {p^q l }, h Y these spheri- 
cal harmonics and the intertwining operator defined by this representation 
is the same as for the original representation. These spherical harmonics 
can be pulled back from the unit X-sphere to the considered X-sphere by 
the central projection ir : X — > X u = X/\X\. Then, for any function 
c(V), function c(Y)tt* (a Vu { Pi qi y) defined on the sphere x ball- type manifold 
is the restriction exactly of those functions (j){ puq .}(\X\, V)o~v u { Pi , qi } f° r which 
4>r p ^ q .y(Rx , V) = c{V) holds. Thus we have 

Theorem 4.6. Both in the fixed and changing basis cases, intertwining 
operator kq induces a well defined action on functions defined by restrictions 
from the ambient manifolds onto the spherex ball submanifolds. This induced 
operator, kq, is well defined for all versions and is the same as the operator 
constructed by a basis of the space of spherical harmonics. 

This induced operator intertwines the Laplacians defined on the spherex - 
ball-type submanifolds. 

The proof of well-definedness for a basis of the space of spherical har- 
monics is the same as on the ambient manifold. The Laplacian on the 
submanifold differs from (|3|) just by the terms Ax and \X\ which should be 
exchanged for A$ x (which is the Laplacian on the X-sphere) and Rx (which 
is the radius of the X-sphere), respectively. Note that the X-directional 
derivatives included into M concern directions tangent to the sphere. Thus 
this term is the same as for the ambient space. In other words, in order to 
have the Laplacian on the submanifold, just the radial Laplacian, A r , of the 
X-space should be dropped from ([3]). Since both the radial and the com- 
plete Laplacians are invariant under the action of the ambient intertwining 
operator, also the Laplacian on the submanifold is invariant. 

5 Domain and range of kq. 

The domain and range of the intertwining operators is determined by a 
function transformation which is noteworthy also without this application. 
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5.1 The dual Radon transform. 



This transform was first investigated in [Sz6j, pages 264-266, where it is 
called boomerang transform. The results provided there include also an 
inversion formula, which, by a new proof, was reestablished by A. Kurusa 
[KuH IKu2| . He called the operator itself dual Radon transform which name 
better describes the area this transform belongs to. We adopt this name, 
however, the following review proceeds with the author's original ideas. 

Let go(r) be a half-line parameterized by arc-length r which has its end- 
point, corresponding to r = 0, at the origin O of M 1 and which is pointing 
to the point 9 6 Sq~ 1 {1) of the unit sphere S'q _1 (1) C M z around the origin 
O. Then (9,r) serve as polar coordinates for the points, Z, of K. 1 . These 
denotations indicate that this transform will be used on the Z-space of El- 
type groups. If f(9, r) is a continuous function defined on M 1 , then, for each 
fixed #o> it determines a cylindrical function fg {Z) defined on the unique 
half-space whose perpendicular projection onto the line spanned by ge ( r ) 1S 
equal to this half-line. If the projection of Z is the point having the polar co- 
ordinates (6*o, r), then, by definition, fg {Z) = f(0o,r). By considering this 
construction for each 9, one can define the function- valued function 9 — ► fg . 
The dual Radon transform, / — * f T , is defined by the integral 

fr ■= [ fed9, (24) 

which can be written also in the form 

fr(Z) := / f(9,(9,Z))d9. (25) 
J(e,z)>o 

Apply Thales' theorem to the last formula to see that the transform is 
defined by the integral of / on the sphere of diameter [0, Z\ by the measure 
d9. Note that this measure differs from the canonical measure of the Thales 
sphere. By (I24p and Fubini's theorem we have: 



Lemma 5.1. Let f(Z) be an arbitrary continuous function and let fi be a 
continuous function with compact support in M. 1 . Then the integral formula 

/ f T (Z)n(Z)dZ = / / f{9,r) m {9,r)drd9 (26) 
Jw Js 1 - 1 Jo 

holds, where hr{9,v) is the Radon transform of fi defined by the integrals 
of this function on the hyperspaces intersecting gg at the points having the 
polar coordinates (9, r) perpendicularly. 
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Formula (|26|) reveals that the considered transform is dual to the Radon 
transform, indeed. The main result in this section is: 

Theorem 5.2. Let f T {d,r) be an arbitrary function of class C 2m with com- 
pact support in M. 1 , where I = 2m + 1 is odd or I = 2m is even. Then f T has 
an inverse, f , regarding the dual Radon transform, which is of the form 

f = J^((frh) {2m \ if l = 2m + l, (27) 

where (2m) means the 2mth derivative of the functions with respect to r, 
resp. 

/= ( " 1 (2^ 1)! ((/ - )fl)[2m ' lf ^ = 2m ' (28) 
where ip[ 2m ^(r) is defined for a function ip(t) on K by 

f°° 1 

V [2m] (r):= ^{<p(r + t)+ip(r-t) (29) 
J o 1 

f 2 f 2m-2 

-2[p{r) + L^(r) + • • • + ^_^^)(r)])(ft. 

Proof. Let /i be a function of class (7 2m with compact support in R 2m+1 . 
From ([26]) we get 



/ 

JM. 1 



{((f T )R) ,2,n, ] T (ZMZ)dZ = (30) 



/ 



C_i "i" 1 

, )^((/.)«) (2m) (^r)^(^r)dr^ = (31) 
gi-i 7 (2vr)^ m 

I (fr) R (0,r)(}-^(^ m \e,r))drd9= (32) 
s 1 - 1 Jo { Z7r ) 

( -i\m p 

f T (z)(^^(^ m) ) r (z)dz = f T (zMz)dz. 



In the last step the well known Radon inverse formula is used. Since \i is 
arbitrary, formula (|27p is established. Formula (|28|) can be established in 
the same way. □ 

Let it be mentioned that all non-trivial isospectrality examples con- 
structed in this paper arise from odd dimensional Z-spaces. Thus, only 
formula (l2"Tl) applies to these cases. 
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5.2 The domain of the intertwining operators. 

The above theorem is used to prove that the function space generated by 
functions of the form J r Q^ p . q .j((f))(X, Z) contains all functions P(X)f(Z), 
where P(X) is a complex valued polynomial and f(Z) is a smooth function 
of compact support on the Z-space. Thus, by using appropriate limiting 
procedures, the whole standard L^-Hilbert space on the (X, Z)-space can 
be generated in this way. Note that the Z-Fourier transform of <p is "twisted" 
with the polynomials appearing in the formula which depend, beside X, also 
on V u . By this reason, it is called also twisted Z-Fourier transform. It is 
this feature what makes the constructions of the above functions highly 
non-trivial. 

The proof of the above statement needs some preparations. For a fixed 
unit Z- vector V® and positive number 5, let T$(V®) be the tube of radius 
5 around the half-line g v o(r). By the standard definition, it is the union 

of those discs, Dg l 2 \r), of radius 5 about the points of the half-line which 
intersect the half-line perpendicularly. The characteristic functions of this 
tube and the half-line, defined on the whole Z-space resp. line determined 
by g v o(r), are denoted by Xsv°(Z) and Xg v o(t) respectively. Then, 

k 12 

J™ ^-^^f^}(xsv S 4>)(x, z) = J] zg.(x)z^.(jr)L^(^)(Z), 

where function L Vu (cf) V o)(t) = X^oW^ivJXs^^ )^). defined on the 

whole line spanned by V® and parameterized by t satisfying t(V u ) = 1, 

vanishes for t < and it is the Laplace transform of <f) v o (r) defined on the 

half- line g v o(r). In the latter formula, this function is described in terms 

of the 1-dimensional Fourier transform Fou+vo defined on the whole line 
n 

spanned by V£. 

Note that the above function does appear as a product of X- and Z- 
depending functions. In the next averaging process they are used to con- 
struct P(X)f(Z) such that one considers the same polynomial P(X) for 
all V u and, in the end, the / appears as the dual Radon transform of an 
appropriate function defined on the Z-space. 

For a complex valued polynomial, P(X), let Pyo(X) be a representa- 
tion of the polynomial in terms of the complex structure Jyo. Further- 
more, let Pyo(X,V u ) be the function defined by replacing V® with gen- 
eral V u in Pyo{X). Therefore, Pyo(X) = P v o(X,vf) holds, but for other 
VuS there are other polynomials defined. Denotation JFqp means that 
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[L=i z vo-(X)zvo-(X) is replaced by P v o(X, V u ) in the above formulas. Then 
we have: 

i im n / tW^V (XsvsM*, Z ) dV u = (33) 

S^oJ s i-i Vol(D ( s ') u 

P(X) [ L c v0 (cf> v o)(Z)dV° = P(X)(h ±v M))r(Z), 

where l^±y u (4>) denotes the function defined on the whole line spanned by 
±V U by the Laplace transforms of functions 4>v u { r ) resp. cf)-v u (r). Thus we 
have: 

Theorem 5.3. For given polynomial P(X) and smooth function f(Z) of 
compact support the product P{X)f(Z) is limit of convergent sequences of 
functions belonging to the domain of an intertwining operator. This sequence 
is constructed by the above method, where function <f>(V) = (/ T - i)(V) 
is derived from (f> by the inverse dual-Radon resp. 1-dimensional Laplace 
transforms defined above on the corresponding half-lines. (The inverse for- 
mula for the Laplace transform is called Mellin's formula. An alternative 
version is the so called Post's formula.) 

The same proof yield those versions of the theorem when function f is 
of the form f(\X\,Z) such that, for any fixed \Xq\, function f(\Xo\,Z) is 
of compact support on the Z-space, or, when this problem is considered on 
a spherex ball-type domain and P(X) is replaced by its restriction, P(X), 
onto the sphere and f(Z) is the same function as before. In these cases 
P(X)f(\X\,Z) resp. P(X)f(Z) are in the domain of the intertwining op- 
erator. In both cases, functions (f>(\X\,V) resp. <fi(V) can be found by the 
same inverse operations defined on the Z-space. 

6 Intertwining of the boundary conditions. 

The most important tool applied in establishing the intertwining of the 
boundary conditions is a theory developed for one- and two-pole functions. 

6.1 Formulas for one- and two-pole functions. 

For a unit Z-vector Z$ and Q S denotation Xqz means that this 
X- vector is in the subspace spanned by Q and Jz (Q)- On the plane 
P{Qi Jz (Q)) spanned by these two vectors the polar coordinates (\XQz \,a) 
are defined such that a(Q) = 0,a(Jz o (Q)) = tt/2 hold. By the restriction 
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a < ir imposed for all Zq, one has a spherical coordinate system on the 
(I + l)-dimensional space, Sq, spanned by Q and all Jz (Q)- Thus these a 
parameter lines are half circles running in the half-plane P + (Q, Jz (Q)) C 
P(Q, Jz {Q)) bounded by IRQ an d containing Jz (Q)- 

Function Qq can be described by this coordinate system as follows. If 
the orthogonal projection, Xq, of X onto Sq is in P + (Q, Jz (Q)), then 

Q Q (X, V u ) = (Q, X Q ) +i([Q, Xq],V u ) = \X Q \ (cos a + i(Z , V u ) sin a). (34) 

Powering performed in QqQq yield: 

Theorem 6.1. On those vectors, X , whose projections onto Sq fall onto 
a fixed a-half-circle around the origin of the half-plane P + {Q , J z {Q)) , a 
1-pole function J~Qpq{4>) has the form: 

r Qpq {<t>){x,z) = [ <j>(\x\,v)%(x,v u )W Q {x,v u )e< z ^dv = (35) 



p+q 

El* 

s=0 



\ p+q cos p+q - 



asm a 



[ A spq {ZQ,V u ) s cj>{\X\,V)e il > z ^dV = (36) 
Jm 1 



P+Q n 

y2\X Q \P+i co S p + i- s a sin s a(A spq (-i) s d s Zo ) / </>(\X\, V)\V\~ S S Z ^ dV.{2,l) 
s=o ^ 

Function 4> S (X,V) = (p(\X\,V)\V\~ s , whose Z-Fourier transform, 4> s , ap- 
pears as the last integral term of \36\). is derived from <p such that it de- 
pends just on \X\ and the Z-variable. Term behind sin s a is denoted by 
A spq (\X\, Zq, Z). If both \X\ and Z are fixed, then A spq is constant for those 
X 's which project onto the half circle determined for Zq by the parameter- 
range a < 7r. On Sq, whose points are denoted by Xq, this function appears 
in the form 

p+q 

\X Q \ p+q co S p + i- s a sin* a A spq (\X Q \, Z , Z), (38) 

s=0 

where, for fixed values of \X\ and Z, the A spq is an s th -order polynomial 
which can be described in terms of the unit vectors Jz (Qu), where Q u = 
Q/\Q\, as follows. 
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Originally, this polynomial can explicitly be determined on the Z-space 
by the expansion {Zq,V u ) s = Ylj=o Bj&Zo' 1 (V u ) in terms of the spherical 
harmonics a s ^ 3 (V u ). Since, for any fixed Zq, function {Zq,V u ) s defined 
on the unit Z- sphere is radial about the center Zq, thus also the spherical 
harmonics are radial about Zq and the convolutions with them are nothing 
but the projections onto the corresponding subspaces of spherical harmonics. 
Thus, 

A spq (\X\,Z ,Z) = J2B j ^ s -ti(\X\,Z ,Z), 
j 

where ip( B ~fi is the corresponding spherical harmonics appearing in the ex- 
pansion of (p(\X\, Z,V U ) = J °° A spq (j) s (\X\, V u , r)e ri ( z > Vu 'dr which function 
is defined, for fixed Z and \X\, by integrals with respect to dr defined for the 
polar coordinate system (V u ,r). 

But this function depends on X. In its final form, it can be viewed such 
that the function determined on the Z-space defines, first, a O-homogeneous 
function on the equator plane Eq u spanned by the X-vectors Jz (Qu)- Then, 
it extends onto Sq such that, on an Xq, it takes the value determined by Zq 
if and only if Xq G P + {Q, Jz (Q))- 

For other points, which are outside of Sq, the function is determined by 
this function and projection onto Sq. Note that the A spq is defined for X 
and not for Xq, meaning that, instead of \Xq\, function (ft involves \X\ to 
this term. The latter parameters stand in front of the formula and are in 
connection with the trigonometric polynomials. 

Such formulas can be established also for TiJ r Q pq ((ft). Functions Qq resp. 
0g are homogeneous harmonic polynomials of the X-variable, thus in cases 
satisfying p = or q = 0, the function in (f35l) is nothing but JiT Q pq ((ft) ■ If 
pq ^ 0, there are new terms, |X| 2r Og r Bg r , appearing in the X-harmonic 
polynomial ILx(&qG>q)- For each r, an additional sum shows up both in ([55]) 
and (|37p . Comparing the first and the r th sums, the p+q is exchanged for p+ 
q — 2r, which is the greatest possible value for s r in the sum. Such a new sum 
can be combined with the first one, where r = 0, by multiplying the r th sum 
by 1 = (cos 2 a+sin 2 a) r . Thus one gets trigonometric polynomials appearing 
in the first sum. By collecting the terms belonging to the same trigonometric 
polynomial, the first term behind the integral sign of (J36j is exchanged for 
a polynomial of the form P spq (\X\, (Z , V u }) = X^=o ^fc r) |^r(^o, V u ) s ~ r , 
resulting the integral terms 
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Pspq(\X\, Zq, Z) 




p spq (\x\,z , v u )4>{\x\,v)S z ^dv, 



(39) 



s 



P spq (\X\,d s z - r J s _ r ) = 2>fc r) (-i) a-r |*r0£ , >--r 



r=0 



behind sin s a in formulas (j36H resp. (I37p . Note that constants A. 
built up by but not equal to the constants A spq . Thus we have: 



arc 



Theorem 6.2. On Sq, function JiTQ vq ($) appears in the form 



p+q 



ip+q 



COS' 



,p+q-s 



asin s aP spq (\X Q \,Z , Z) 



(40) 



s=0 



where P spq is explicitly described in 139\) . For fixed values of\X\ and Z also 
this term is an X-depending s th -order polynomial which appears in the same 
form as A spq does. But this one has also lower order terms, (Zo,V u } s ~ r , 
beneath the main term. For other points not being on Sq also this function 
is determined by projections onto Sq, in which case it is defined in terms of 
\X\ and not \Xq\. 

The above constructions restricted onto spheres Sr x provide the formu- 
las on sphere x ball- type domains. In this case function \X\ is constant, thus 
functions (f> s - r depend just on Z and Zq. Let it be pointed out again that 
the latter variable is involved by the assumption Xq S P + (Q , J r z (Q)) , i- 
e., it is determined by X over which the Fourier transform in the Z-space 
is performed. In other words, it is an X-depending function whose precise 
denotation would be Zq(X). 

Since the Z-balls, Bji z (X), where X £ Sr x , are naturally identified on 
this trivial ball-bundle, they determine the same functions in the Z-space 
for all those X's which project onto the half-plane P + (Q, Jz (Q))- Thus 
functions !FQ pq {4)) resp. 'HJ : Q P q(4>) appear in the form ([38]) resp. P0|) such 
that \X\ = Rx is constant in this case. 

Later on, we need these functions described also on circles, C = P2^Sr x , 
which are represented as intersections of 2-dimensional linear subspaces, 
P2, with Sr x . If Q G P2, these functions are perfectly described by the 
above formulas also on these circles. Therefore, we suppose Q P2. First 
also suppose that P2 C Sq holds, in which case the computations below 
are carried out on the 3-space, Sq%, spanned by P2 and IRQ. This space 
intersects Sr x at the 2-sphere denoted by Sr x 2- The north-pole, O, of 
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this sphere is cut out by the ray M+Q- The north-pole, Oc, on the circle is 
defined by the closest point to O. Let ate be the angle parameterization of C 
with origin Oc such that on both sides < etc < vr hold. The angle between 
C and the great circle CEo with center O (equator) on the 2-sphere Sr x 2 is 
denoted by ftc- It is uniquely determined by the assumption < (3c — ir/2. 
For a point P 6 C satisfying < oic < vr/2, let Cp be the great circle 
connecting O and P, which intersects CEo at a point N perpendicularly. 
If Mc = C H CEo, then the spherical sine theorem applied to the right 
spherical triangle PNMc yields 

sin(f - a(P)) sin(f - a c (P)) , N , s , , 

- - — cosa(P) = sin /3 C cos a c (P) (41) 



sin Pc sin 



2 



This equation along with sin a(P) = yl — sin 2 Pc cos 2 otc(P) imply that 
functions (|38[) and (I40p restricted onto C can be expressed in terms of 
cos ac{P)- 

It is a very important issue to understand the precise appearance of 
functions A spq and P spq on these circles. They appear as polynomials on 
the equator circle CEo but on C they appear as functions which are pulled 
back from the equator to the C by the central projectivity to ■ C —* CEo- 
This to can be explicitly computed as follows. 

Parameterize both P% and Eo by complex numbers z and z' respectively 
such that these coordinate systems have common imaginary axis M+i C 
P2 U Eo and z = 1 and z' = 1 correspond to P = Oc and iVp = Cp n CEo 
respectively, where Cp is the great circle connecting O and P. This circle 
intersects CEo at ./Vp perpendicularly. Pick up also such unit complex 
numbers, u and u' , between the units and the imaginary numbers which are 
corresponded to each other by the To- 

Actually, the to is a real projectivity between the two projective lines C 
and CEo defined by identifying the antipodal points on the great circles C 
and CEo- Thus, it can be described in terms of the real cross ratio defined on 
these projective lines. But this real one is the same as the complex cross ratio 
defined on the complex planes if the points are laying on the same half-circle. 
(This statement is well known in conform geometry of 2-spheres.) Also note 
that the common imaginary numbers cut the great circles into half-circles 
which are corresponded to each other by the to, therefore, this projectivity 
can be described in terms of the complex cross ratio by the relation (z' = 
to{z), 1', i', u') = (z,l,i,u). This equation describes the to as a fractional 
linear function (Mobius transform) of the form z'(z) = (az + b)/(cz + d). 
Since such a function preserves the circles and the three corresponding points 
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(l',i',u') and (1, i, u) are on CEo and C, the transformation is really denned 
between the two circles. Thus we have: 

Lemma 6.3. The to is a fractional linear function (Mobius transform) be- 
tween the two great circle which pulls back a trigonometric function cos u a = 
(l/2)(uz + uz), defined on the complex plane Eq by a fixed complex unitu, to 
the trigonometric rational function (acos u +b)/(ccos u +d) defined on C . All 
trigonometric polynomials can be generated on CEq by the functions cos u , 
therefore, all trigonometric polynomials are pulled back to a trigonometric 
rational function defined on the other great circle. Functions A spq and P spq 
can be described by the trigonometric functions ip( s ~^ defined on the unit 
sphere of Eq. Thus they are trigonometric polynomials on CEq which pull 
back to trigonometric rational functions defined on C . 

For a P2 which is not subspace of Sq these functions can be determined 
by projecting it into Sq. Almost every plane projects to a plane of Sq. If P' 2 
is the projected plane, then the sought functions on P2 are the pull-back's 
of the corresponding functions defined for P' 2 . Note that cosa^ arises from 
a linear function, therefore, so does the pull-back function whose kernel is 
the pull-back of the line (linear subspace) connecting the points Mqi and 
—Mqi defined above for C . Then the closer midpoint, Oc, to the O, which 
is between the pull-back-points Mq and —Mq on C, is called the north- 
pole on C. This point determines the parameterization ac- The functions 
restricted onto C are described in terms of this parameter. Thus we have 

Theorem 6.4. When the constructions are restricted onto a fixed sphere 
Sr x , functions TQ vq {$) resp. riTQ pq {^>) appear in the form $38\) resp. (Jty 
such that functions A spq and P spq involve the constant \X\ = Rx- In the 
following statement denotation R spq can be replaced by any of these two 
functions. 

On a circle, C = P2CiSr x , represented by intersection of a 2-dimensional 
linear subspace P2 with Sq, these functions are of the form 

p+q 

Y,K p c +q (smp c > cosa c ) p+q - s (l - sin 2 C > cos 2 a c )^R spq (Z , Z), (42) 

s=0 

where C is the projected circle cut out by the projected 2-space P 2 which 
intersects the equator E' Q on the projected space Sq^ at angle f3c> ■ Constant 
K^ q is due to the fact that the pulled back linear functions are restricted to 
a circle in this process. For circles in Sq this constant is Kq = Rx- 

For any fixed Z and \X\, function R spq (Zo, Z) defines a trigonometric 
rational function on C which is the pull back of an s th -order polynomial 
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with such a combined map, where the first map, t'q, is a Mobius transform 
between C to CE'q and the second one takes C onto C by a projection. 

We need these theorems in the more general case when, instead of Q, one 
considers a pair, (Q (a) ,Q ib) ), of vectors and functions <p(\X\,V), Qq Vu ®qv u > 

are exchanged for the following ones 

respectively. Such functions are called 2-pole functions which can be inves- 
tigated in two ways. They can be considered either on subsets (Xp \ X^) 

defined by a fixed X^\ or, on the similarly defined subsets (X^ a \ Xp^). 
Because of the exact similarities, only the first case should be described, 
when, functions 

<Hi4 a) i,ix( 6 )|,n ®gw»®QWn' n K ?o) (") ( 43 ) 

depend (non-trivially) just on \X^\ and V. On a circle (X^\ C^) the con- 
sidered functions appear in the form (I42D where the last function is defined 
by those listed in (|43p and the other functions are defined on \~( b \ 

6.2 Intertwining of the Dirichlet conditions. 

The Dirichlet Intertwining Theorem will be established, first, for a constant 
basis, Qf. The changing basis case will be traced back to this first one. 

Observe that functions cos p+q ~ s (a) sin s (a), satisfying < s < p + q are 
linearly independent, furthermore, for any fixed Zq and Z, function A spq is 
constant on the a-parameter line determined by Zq. These two statements 
yield the following theorem obviously. 

Theorem 6.5. A function J-Q pq {4>) satisfies the Dirichlet condition at the 
boundary points (X, Z) if and only if functions A spq vanish on the sphere 
Sr z , for all Zq(X) and < s < p + q. Regarding HT Q P q((j>) , this condition 
is P spq = 0, for allO < s < p + q and Zq(X) at any boundary point Z £ Sr z . 

For fixed Q and natural numbers p and q, function spaces &Q pq resp. 
^Qpq are defined by the L 2 function spaces spanned by functions of the 
form !FQ pq {4)){X, Z) resp. TLJ- 'Q pq (cp)(X , Z) , where cf>{\X\,V) (which de- 
pends, non-trivially, just on V on spheres sphere-type manifolds) can be an 
arbitrary L? -function. For fixed Q but running p and q, all these spaces 
sum up to the total space *q = Y, p , q *Qp? = H <2 = J2 p , q s Qpq- Then, for 
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functions ipQ resp. ip'g = kq f (pq) from <J?q resp. <frg the Dirichlet condi- 
tion is satisfied always simultaneously. Actually, the intertwining operator 
kq f : <&q — ► <&g between these total spaces is induced by a point transforma- 
tion of the form (TQ(X),idz), where the Tq is an orthogonal transformation 
on the X-space, depending on Q. 

Such simple proof can be given only for total spaces denned by a fixed 
pole Q. The proof is much more difficult on the complete L 2 -Hilbert space, 
which can be represented both as E Q es P an R (Q F ) *Q> and EQes P an R (Q F ) *'q. 
i. e., by the sums of all functions defined by all poles, Q, which are in the 
real span of Qp. The main idea of such an extension is as follows. 

Suppose that a function (p(X, Z) = ^2q. pq^X, Z) satisfies the Dirich- 
let condition. Decompose each Qq in the form @q(X,V u ) = (Q,X) + 
i(J Vu (Q),X) &nde Q (X,V u )(Q,X) - i(J Vu (Q), X), which, after multiplica- 
tions, result the decomposition 

<P = Pevnj + foddj = ¥Q l evn J + VQioddj > (44) 

Qi Qi 

where the first function involves all terms having even number of Jy u while 
for the other one this number is odd. For a fixed boundary point Z, these 
functions depend just on X. One can prove, by formula (|42p . that these two 
functions are in two completely independent subspaces of functions. The 
proof will be based on the fact that, on a circle C, the first function appears 
as a trigonometric rational function depending on sin 2 0c cos 2 ac while the 
other function depends irrationally on these terms. This independence im- 
plies then that both p evnj (X, Z) and (p a ddj (X, Z) must satisfy the Dirichlet 
condition. In the next step one decomposes these functions in the form 

— '-P evn, j evn j ^ (X, Z) + p> evn j odd ^ 

(X,Z), (45) 

Poddj(X, Z) 

V odd j evn j 

(X,Z), (46) 

where the first term of each function involves all terms having even number 
of Jy \ while, for the second one, this number is odd. Now using the double 
pole version of (|42p . one can see that all 4 functions in this final decom- 
position fall in completely independent subspaces. But, then, all these 4 
functions must satisfy the Dirichlet condition. Now we observe that 



^Qf ■ ¥ 'evn j evn ^ * 'P evn j evn j ^ 

(X,Z), (47) 

k Qf • ¥ 'evn j odd 

(X,Z), (48) 

k Qf : Poddjevn ^ *■ foddjevn ^ 

(X,Z), (49) 

K Q F '■ Poddjodd ( b) < -Poddjodd J ( b ) 

(X,Z), (50) 
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which relations are due to = = - j( ft ). Thus <p and ^ satisfies 

the Dirichlet condition simultaneously. 

For completing this proof only the above mentioned Independence The- 
orems should be established. Suppose the contrary, there is a function (p 
which can be represented as linear combinations both of evnj- and oddj- 
type functions: 

S ^Qrevnj =f = Yj fQioddj ■ ( 51 ) 

Qr Qi 

We may suppose that the trigonometric polynomials cos Pj+qj ~ s sin s are of 
the same order n = p j + qj and there are only finite linear combinations in 
this expression. Consider a circle C on which each term in the linear combi- 
nations appears in the form (|42p with the corresponding constant sin /3qi and 
origin (north pole) Ocj on C. (The origin (pole), Oc, on C is constructed 
above Theorem 16.41 ) 

Next we work just on the right side of (|5ip. involving only the oddj- 
type functions. They can be sorted out into classes according their pa- 
rameters sin P c > and Ocj- Re- numerate these classes with m in the form 
{(f3c>m, Oc q )\m = 1,2, ... ,d} such that at least one of the parameters is 
different for two distinct m's. From the partial sum determined by a class 
factor out (1 — sin 2 Pc'm cos 2 acm) 3 • Thus, it appears in the one-term-form 
S m (l — sin 2 flc'm cos 2 »Cm) 5 i where S m is a rational trigonometric function. 
Furthermore, 

d d,oo 

<P = Sm\J 1 - sin 2 Pc'm cos 2 a Cm = ^2 A r S m sin 2r f3 C ' m cos 2r a Cm , 

m=l m,r=l 

(52) 

where YlrLi A r x r is the Taylor expansion of y/1 — x about x = 0. All the 
coefficients A r are non-zeros. Since the left side of (|5ip is a trigonometric 
rational functions, all terms 

d 

^ ^ Sin 2r /? C / m COS 21 " QCm, (53) 
m=l 

considered for a fixed r, vanish for big numbers r, say, if r > N . 

Now observe that there exist an everywhere dense open subset, B C C , 
such that, for any fixed point P £ B the function values sin Pc'm cos ac m (P), 
considered for all m = l,...,d are distinct and, therefore, the d x d- 
matrix (T mk = sm 2( ~ k+N ^ {3 C 'm cos 2( - k+N ^ ac m {P), where k = 1, . . . ,d, is non- 
degenerated. Thus equations ^2 m T m kS m (P) = 0, V/c, implies S m (P) = on 



31 



the whole circle C. This argument can be repeated in those cases when S m 
is substituted by S m sin 2r flc'm cos 2r ctcm , Vr > 0. Therefore, ip = 0, which 
proves the independence of the considered function spaces completely. 

This argument repeated for 2-pole functions proves the desired inde- 
pendence, first, on the subsets (X^ \ X^). But this statement obviously 
implies the independence on the whole X-space. 

Both independence theorems can be formulated in terms of polynomials 
such that, after considering the decompositions, the evn j- resp. odd j-type 
subspaces are spanned by functions involving even resp. odd number of J's. 
This observation leads to a simple establishment of the independence and 
Dirichlet-intertwining theorems regarding changing basis cases as follows. 
First consider the complex matrix Cij(V u ) transforming the fixed basis Qf 
to the changing one, Q(V U ). It is obvious that, both in the 1-pole and 2-pole 
cases, both type of subspaces regarding the two systems are transformed to 
each other by the non-degenerated map 

u '■ jP QHp;,<ft} ~* ^QiVuHpuqi} ' n:F Q F { Pl , qi } MFQ(V u ){ Pi ,qi} (54) 

induced by the basis-transformation. This proves both independence theo- 
rems for the changing basis case immediately. Thus we have: 

Theorem 6.6. The kq intertwines the Dirichlet condition both in the fixed, 
Qf, and the changing basis, Q(V U ), cases. 

The proof is based on the Independence Theorem stating that the 
total space 

^Q F ,n = ^2 ^QfIr,?;}' 
{n=J2(pi+q l )} 

defined for a fixed n and Z , is the direct sum of the independent subspaces 
^Q F ,n,evnj and J~q f t n,oddj ? where, after implementing the above described 
natural decomposition, the functions from the first resp. second space con- 
tain even resp. odd number of J 's. Both of these subspaces further decom- 
pose into the independent subspaces Fq f , n ,par j ,evn j{b) and Fct F , n ,parj,odd j{b) 
defined by the options par = evn or odd given for the parities of the number 
of J( fc ) 's in the expressions. Vector Z should not be the same for the par- 
ticipating functions but it can be chosen individually and independently both 
for the even- and odd-type functions. 

The independence guaranties that all 4 component functions of a ip and 
ip itself satisfy the Dirichlet condition always simultaneously. The same 
statements hold for the function spaces TtTQp^n o,s well as for both versions 
of function spaces defined by changing basis fields. 
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7 Intertwining of the Neumann conditions. 



The Neumann conditions create a new more complicated situation which 
requires reformulations of the proofs given for the Dirichlet conditions at 
several points. The Z-Neumann conditions, however, which require the van- 
ishing of the derivatives of functions taken from the Z-radial directions at the 
boundary points, can be strait-forwardly traced back to the Dirichlet condi- 
tions. Let it also be mentioned that the proofs on the boundary manifolds 
exploit only the intertwining of the Dirichlet and Z-Neumann conditions. 
By this reason, the Z-Neumann conditions are considered first. 

In the following computations the integral defining the Z-Fourier trans- 
form is considered on the polar coordinate system. The computations are 
carried out by formulas d\ z \e l{Z ' V) = i{Z u , V U )\V\S Z ^ = \Z\- x \V\d\ V \£W 
combined with integration by parts. Without loosing the generality, one can 
suppose that the test-function, <f>, vanishes at the infinity. Then, in terms 
of 4>' := d\z\<j>, we have 

d\z\FQ{ Pm M = -\Z\{^Q{ Pi ^\ V W)+ l ^{PmM)- ( 55 ) 

Therefore, a function •^ r Q{ Pi ,g i }('/ ) ) satisfies the Z-Neumann condition if an 
only if FQ^p.^ydVlcj)') + lJ r Q{ Puq .}((j)) satisfies the Dirichlet condition. Since 
the Dirichlet condition is intertwined in all cases, we have 

Theorem 7.1. The Z-Neumann condition is intertwined both in the fixed 
and changing basis cases. 

From now on, the standard Neumann condition is scrutinized. By for- 
mulas (3.2) and (3.7) of |Sz2j . the normal vector at a boundary point (X, Z) 
and the Laplacian on the boundary manifolds are of the form 

fi = A(\X\,\Z\)X U + B(\X\,\Z\)J Z (X) + C(\X\,\Z\)Z U , (56) 

1 l ~ l 
A = A Sx{z) + (1 + -\X\ 2 )A Sz{x) + £ d a D a ; (57) 

a=l 

where Sx(Z) is the X-sphere over Z and Sz{X) is the Z-sphere over X, 
furthermore, {di, . . . is an orthonormal basis in the tangent space of 

the Z-sphere Sz(X) at Z. On sphere x ball- type manifolds the first term of 
([56]) should be omitted. 

The standard Neumann condition is considered, first, for one pole func- 
tions. The derivative with respect to the normal direction is built up by 
X- resp. Z-radial-derivatives and Jz(X)». If the foot of perpendicular 
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through X to Sq is Xq, and thus X = Xq + Xq holds, then J Z {X)» = 
Jz(X Q ). + J z (X^): First J z (X Q )m is considered. If X Q = \X Q \(cos(f3)Q + 
sm((3) J Zo (Q) e P(Q,Z ), then 

J Z (X Q ) = \X Q \(cos(f3)J z (Q)-sm(P)(Z ,Z)Q + sm(P)J z ±J Zo (Q)), (58) 

where Z -1 - is the perpendicular component of Z to Zq. Therefore, the J z ±. 
and J Zo are anti-commuting and the J z ± J Zo is a skew endomorphism. Thus, 

(Q, Jz(X Q )) = -\X Q \ smp(Z , Z) , (59) 

(J Zo (Q),J Z (X Q )) = \X Q \ cos p(Z ,Z), (60) 

Jz(Xq) • cosa = -\X Q \ smf3(Z , Z) , (61) 

J z (X Q ).sma = \X Q \ cos (3{Z ,Z), (62) 

Jz(X Q ) . HrQ„{<l>){X, Z)=J2 \X Q \ p+q cos"^"'^) sin s (/?)Sj OT , (63) 
w/iere Sj* = -|Xq|(Z , Z)((p + g - s + l)P (s _ 1)p , - (s + l)P (s+1)w ). (64) 

The computations with J z (Xq)» are based on 

{Q,J z (X^))=0 , {J Zo (Q),Jz(Xq)) = -{Jz±Jz (Q),X%). (65) 

Since differentiation J z {Xq)» acts, non-trivially, on the considered functions 
only by its contribution (J Zo (Q), J z {Xq)){cos (3d a + sina5 r ) to the d a - and 
the radial <9 r -direction, therefore: 

J z (Xq) »cosa = (J z ±Jz (Q),Xq) cos ft sin /?, (66) 
Jz(X^) . sina = -(J z ±Jz (Q),X^) cos (3 cos (3, (67) 
Jz(X^).HF Qpq (<P)(X,Z) = |X |^^cos^- s (/3)sin s (/3)5^, (68) 

s 

S s - L p? = -(J z xJz (Q),^)((p + ? )|X Q |- 1 sin/3P s?)? -cos/3^), (69) 

-Dspg = (p + q - S + l)P(s-l)pq ~ (« + l)P( a+ i)p,. (70) 

A preliminary version of the standard Neumann condition for a one- 
pole function can be stated in the following form. A ^-generated one-pole 
function satisfies the Neumann condition if and only if 

R spq = (Ad ]x] + Cd\z\)P spq + B{S T spq + ~Si pq ) = (71) 

holds at the boundary points, for all Zq{X) and < s < p + q, where co- 
efficients A, B and C are defined by the normal vector fi. Already this 
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version reveals that only the intertwining regarding the first two terms 
in this condition can be traced back to the intertwining of the Dirichlet 
conditions. The other two terms involve functions such as {Zq,Z) and 
(Jz (Q),Jz(Xq)) = \X^\{J Zo (Q), Jz(Xq )) which appear outside of the 
integral terms. It is noteworthy that the second function is zero on spaces 
H^ a ' b ^ and H^ a,b \ This is due to the fact that the irreducible components 
-ffg 1,0 ^ and H^'°\ yield the well known J 2 -condition, meaning, that for any 
product J z ±Jz there exist J z such that J z ±Jz = J z holds. On arbi- 
trary H-type groups, for fixed Z and unit vector Xq Q , all these terms define 
polynomials which are suitable to establish the independence theorems seen 
for the Dirichlet condition also for the Neumann condition. Similar formu- 
las can be established also for the two-pole functions which also yield the 
corresponding independence theorem. Finally we get: 

Theorem 7.2. The kq intertwines the standard Neumann conditions both 
in the fixed, Qf, and the changing basis, Q(V U ), cases. 
The proof is based on observing that the total space 



defined for fixed Z, n and running {pi,qi} satisfying n = ^2(p% + qi), is a 
direct sum of independent subspaces 'HJ : Q F ^ n ^ vnj and 'HJ r Q F ,n,oddj7 which 
are defined such that the functions from the first resp. second space con- 
tain even resp. odd number of J 's (i. e. , sink's, according to the above 
decomposition) . Both subspaces further decompose into the independent sub- 
spaces fj, • J r Q F ,„, P ar J ,ei,n j(f)) and fx • Fq f , n , P arj ,odd j(b) defined by the options 
par = evn or par = odd, available for the parity of number of 's in the 
expressions. The independence guaranties that all the 4 component func- 
tions of a ip along with ip satisfy the standard Neumann condition always 
simultaneously. The same statements hold for the function spaces Tq f ^ as 
well as for both versions of function spaces defined by changing basis fields. 

8 Intertwining on the boundary manifolds. 

Since the intertwining operator preserves the Dirichlet condition, by re- 
strictions, it induces a well defined bijection between the L 2 spaces defined 
on the boundaries. Since each smooth function on the boundary extends 
to ones satisfying the Z-Neumann condition, furthermore, this condition is 




{Pi,Qi} 
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also preserved by the operator, it is enough to represent the functions on 
the boundary by restrictions of those satisfying the Z-Neumann condition. 

If di = d\z\ is the Z-partial derivative with respect to the normal direc- 
tion Z u , then the angular momentum operator M (resp. M) on the ambient 
(resp. boundary) manifold differ from each other just by d\z\Dz u *- This 
operator vanishes on functions satisfying the Z-Neumann condition, thus M 
and M acting on these functions provide the same results. This argument 
proofs that not just M and M' but also M and M' are intertwined by the 
operator. This is the most crucial part in the proof of the pursued theorem. 
The intertwining regarding A$ x (z) has already been established on the am- 
bient manifold, thus one should consider only Ag z (x) • Since the intertwining 
regarding Az is established on the ambient manifold, only the intertwining 
of the radial Laplacian A\z\ should be established. Since this operator acts 
on functions satisfying the Z-Neumann condition, the question is if d? z , 2 is in- 
variant under the action of the operator. This statement immediately follows 
from the following computations where the integral defining the Z-Fourier 
transform is considered on the polar coordinate system. The computations 
start out with d 2 z? S z ^ = -(Z U ,V U ) 2 \V\ 2 S Z > V ) = \Z\- 2 \V\ 2 d 2 v? S z > v ) 
and are completed by integration by parts. Then, in terms of <j)' = d\z\(j>, we 
have 

dfz^QpM = \Z\-\r Qpq {\V\ V) + 2(1 + l)F Qpg (\V\cf>') + IFqMI (72) 

where the \Z\ = Rz is a constant. Like in case of the Z-Neumann condition, 
this formula establishes the desired intertwining property for d 2 z ^ 2 - 

9 Intertwining on solvable extensions. 

The isospectrality theorems naturally extend to the solvable extensions. The 
Laplacians on the ambient- and boundary-manifolds, furthermore, the nor- 
mal vectors at the boundaries are described in formulas (1.12), (3.30), and 
(3.29) of [Sz2]. The generator functions are of the form (^(|X|,t, V) in this 
case, till, the intertwining operator is defined by the same Z-Fourier trans- 
form like on H-type groups. I. e., the intertwining on the solvable group is 
completely determined by its action induced on the nilpotent group. The 
details in [Sz2] show that the terms due to the t-variable, which are in a cer- 
tain combination with the terms of the Laplacian defined on the nilpotent 
group, make no effect on proving the intertwining for the solvable groups in 
the same way as for the H-type groups. 
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10 The new striking examples. 



There is a subgroup, Sp(a) x Sp(6), of isometries on a Heisenb erg- type group 
H^ a ' b ^ which acts as the identity on the Z-space. Note that these isometries 
act transitively on the X-spheres of H^ a+b ' ^ . 

The complete isotropy group of isometries fixing the origin is (Sp(a) x 
Sp(6)) • 50(3), where the action of SO(3), described in terms of unit quater- 
nions q by a q (X 1 , . . . , X a+b , Z) = {qX{q, ... , qX a+b q, qZq), is transitive on 
the Z-sphere. The elements of this isotropy group induce isometries on the 
sphere x sphere- type submanifolds, furthermore, there is proved in the Ex- 
tension Theorem of [Sz2j that these are the only isometries on these subman- 
ifolds. Note that these isospectral manifolds in a family have non- isomorphic 
isometry groups of different dimensions such that they are homogeneous in 
jj(a+b,o) ^ jj(o,a+b) ^ w j 1 £ e o^gj. me mbers are locally inhomogeneous. 

The sphere x sphere- type submanifolds of the solvable extensions of H- 
type groups is defined such that, over each point of a sphere in the X- 
space, one considers the same geodesic sphere around the origin (0, 1) of the 
hyperbolic (Z, t)-space. The isotropy group of isometries acting on SH^ a ' b ^ , 
where ab / 0, is (Sp(a) x Sp(6)) • SO(3), while it is Sp(a + b) ■ Sp(l) on 
SH^ a+b ' ^ = SH^°' a+b \ The above statements extend also to these groups. 

These are the new striking examples brought by the reconstructed inter- 
twining operator. The sphere-type striking examples discussed in the earlier 
papers can be explained similarly. They are the geodesic spheres defined by 
the same radius for the family SH^ a ' b ^ . In this case the metric on SH^ a+b ' ^ 
is two-point homogeneous, therefore, having homogeneous geodesic spheres. 
The geodesic spheres on the other groups are locally inhomogeneous. 

It should be mentioned that Schiith [Schl] constructed the literature's 
first isospectral metrics defined on simply connected manifolds on Cartesian 
products of spheres. Among them there are also sphere x sphere- type mani- 
folds. Her construction arises from a completely different setting, however, 
and all provided metrics are locally inhomogeneous. 



11 Isospectralities for cr-equivalent metrics. 

A ex-deformation of an endomorphism space J z is defined by an involutive 
orthogonal transformation, a, of the X-space which commutes with all endo- 
morphisms of J z . The er-deformed endomorphism space consists of endomor- 
phisms aJz- By using irreducible decomposition regarding the orthogonal 
Lie algebra generated by the elements of J z , for any cr-deformation, there 
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exist a decomposition v = v( a ) © with components invariant under the 
action of the endomorphisms such that for their restriction onto the com- 
ponents the relations crj^ = Jg and oJ^ = ~J^z > n °ki. Note that the 
family, Jj- a,b \ of Cliffordian endomorphism spaces defined by the same a + b 
and I consists of u-equivalent endomorphism spaces. In papers |Szll ISz2| 
the isospectrality is stated on the ball- and sphere-type domains of such 
a-equivalent 2-step nilpotent Lie groups and their solvable extensions whose 
endomorphism spaces contain at least one anticommutator. Thus the ex- 
tension of the isospectrality theorems to cr-equivalent 2-step nilpotent Lie 
groups and their solvable extensions provide plenty additional examples to 
those produced by the anticommutator technique. For the sake of simplic- 
ity, we consider such two step nilpotent Lie groups whose endomorphism 
spaces contain at least one non-degenerated endomorphism. Then almost 
all endomorphisms must be non-degenerated acting on an even dimensional 
X-space. 

First, we look for the necessary modifications which make the techniques 
developed for H-type groups working also for cr-deformations. The Lapla- 
cians in these general cases differ from the Laplacians of H-type groups just 
by the term (1/4)|X| 2 A Z , what is now (1/4) £(J a (X), J p {X))d 2 ap . Also the 
intertwining operator, defined for the changing basis case, must be modi- 
fied as follows. Let (Qi(V u ), . . . , Qk^iYu)) = (Q > Q ) be an appropriate 
orthonormal basis field such that each Qv u i is an eigenvector of J v with 
eigenvalue —X 2 (V U ). Let Jy u be the normalized endomorphism which has the 
same kernel as Jy u and is defined by (l/Xi(V u )) Jy u on the maximal eigensub- 
spaces belonging to \ > 0. Note that the kernel is trivial for an everywhere 
dense open subset of the unit vectors V u , furthermore, this endomorphism 
may not be in J z . Then, by definition, Qq^ Vu ^(X,V u ) = (Q + i Jy u (Q) , X) . 
The complex coordinate system defined by this basis for non-degenerated 
endomorphisms is denoted by {zy u i, • • • , Zy u k/2\- Then the intertwining is 
defined by these changing complex coordinates such that functions <j> should 
be of the form . . . , |Xy r |, V), where X = Xyi is the decomposi- 

tion regarding the eigenspaces of J v . This requirement is slightly different 
from that what is considered on H-type groups, but one can reach to them 
in the same way: Start with functions depending just on V, first. Then, it 
turns out that the same operator is defined by the above more complicated 
functions. Furthermore, the domain is the largest possible, containing the 
complete L 2 - function space. Also Theorem 14.41 concerning the intertwining 
of the Euclidean Laplacian and radial functions on the X-space, remains 
true for cr-deformations. Therefore, by Jy = aJyoJy = a 2 J v = J v and 
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formulas (|22p . (|23p . which also extend to cr-deformations, this is indeed an 
operator intertwining the Laplacians A and A' term by term. 

One should check out also the intertwining of the boundary conditions. 
First note that the technique developed for H-type groups works out strait- 
forwardly only on groups where the intertwining can be established also 
with a fixed basis Qf, therefore, one can use one- and two-pole functions 
with poles being in span^Qp). Also well-defined polar coordinate systems 
established on IRQ © J Z (Q), where J z is spanned by endomorphisms of the 
form Jy u , are necessary conditions for this technique. All these requirements 
are satisfied only in those cases where, for all unit pole Q G spariM.(QF), the 
unit vectors J Z (Q) form an everywhere dense open subset of the unit sphere 
of the space spanned by these vectors (equator). Then, any such vector is 
connected with the pole by an a-parameter circle defined for < a < ir such 
that it has the parameter ir/2. The possible missing circles, which are due 
to the degenerated endomorphisms, can be implemented by limiting. The 
extension to this cases works out after other additional modifications. 

First we check on formula (|34p of Section 16.11 where, on a parameter- 
circle, the corresponding Zq should be exchanged for Z*(Zq) defined by 
the dual of the functional tp(Z) = (Jz (Q), Jz{Q)) ■ Yet, the Z*{Zq) is 
a polynomial function of Zq, implying that functions R spq , introduced in 
(|42p . will be polynomials on the above equator. Thus the computations 
can be processed in the same way as earlier. Since for cr-deformations the 
relations <p(Z) = (p'(Z) , Z* = Z^ hold, the proof regarding the Dirichlet 
and Z-Neumann condition can be completed by the same argument seen for 
H-type groups. 

Regarding the Neumann condition the Zq inside of the integral term 
should be exchanged for Z*, while terms (Zq,Z) resp. (Jz^JzoiQ), Xq) 
outside of the integral should be exchanged for much more complicated ex- 
pressions. Even so, they provide polynomial functions and the modified 
computations concerning formulas ([59]) - ([64"]) along with (Jz*Jz (Q),Q*) = 
{J' Z *J Z (Q),Q*) yield the intertwining also of the Neumann conditions for 
the cr-deformations whose endomorphism spaces satisfy the above condi- 
tions. In these cases, the theorem extends also to the boundary manifolds 
and the solvable extensions. 

Fortunately, one should not go through the steps of this complicated 
construction which is incomplete without scrutinizing the question of the 
existence. Basically, what the above proof exploits is that the operator 
introduced by the changing-basis-technique for cr-deformations intertwines 
both the Laplacians and boundary conditions if the manifold satisfies the 
independence theorems. But this theorems are certainly yielded on groups 
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having endomorphism spaces which are the results of "small" perturba- 
tions performed on the endomorphism space of H-type groups such that, 
by choosing a new endomorphism space close to the Clifordians, one per- 
turbs the endomorphism space acting on the irreducible space M. r ( l \ This 
defines uniquely determined perturbations for the reducible endomorphisms 
(see more details in the end of Introduction). If the perturbation is "small" 
which changes the endomorphisms just slightly, the subspaces in the inde- 
pendence theorems (which are closed in the ambient Hilbert space) keep 
being independent. Even the conditions for the existence of polar coordi- 
nate systems are satisfied on manifolds defined by "smaller" perturbations. 
However, the above construction can be left out completely because the idea 
of perturbation provides more examples than those provided by the above 
process. Also the non-isometry proofs, guaranteeing that the considered 
isospectral metrics have different local geometries, are inherited from the 
metrics being "slightly" perturbed. For the latter metrics the non-isometry 
proofs are completely established in [Szl, Sz2]. Thus we have: 

Theorem 11.1. The intertwining operator introduced for a -equivalent groups 
by the changing-basis-technique always intertwines the Laplacians both on 
ballx ball- and spherex ball-type manifolds. Furthermore, there exist an open 
neighborhood, U , in the space of l-dimensional endomorphisms spaces with 
endomorphisms acting on the irreducible space W^ 1 ' of a given Clifford endo- 
morphism space, J z such that the latter one is in U and all 2-step nilpotent 
groups and their solvable extensions constructed by endomorphism spaces be- 
longing to U satisfy the independence theorems. Therefore, for these groups, 
the intertwining operator intertwines also the boundary conditions. In other 
words, small perturbations performed on the endomorphism spaces of H-type 
groups provide a wide range of a -equivalent groups which are isospectral on 
the corresponding ballx ball- and spherex ball-type submanifolds. 

This theorem extends to solvable groups as well as to the boundary man- 
ifolds, both in the nilpotent and solvable cases. For dimensions satisfying 
I = 4ra + 3 also such U exist for which the corresponding metrics in an 
isospectrality family have different local geometries. 

Remark. The key point in the above process dealing with one-pole func- 
tions is that the multi-linear function 

H(X,X*,Z,Z*) := {J Z (X), J Z *{X*)) (73) 

does not change during cr-deformations. The Ricci tensor can be described 
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in terms of this function by 

R(X,X*) = -(l/2)Y J H(X,X*,Z a ,Z a ), (74) 
R(Z,Z*) = (l/A)Y,H(E i ,E i ,Z,Z*) , R(X, Z) = (75) 

(cf. formula (1.9) of |Szlj ). thus also this tensor is not changing during 
cr-deformations . 

The Gordon- Wilson [GWJ isospectrality examples were constructed on 
ball x torus- type manifolds by spectrally equivalent endomorphism spaces, 
meaning the existence of orthogonal transformations associating isospectral 
endomorphisms to each other. More precisely, they constructed continuous 
families of metrics which are isospectral on functions. These metrics are not 
isospectral on 1-forms, however, due to the fact that the norm of the Ricci 
tensor is changing during these deformations |^ch2j. The question arises if 
the domains investigated in this paper are isospectral on the Gordon- Wilson 
examples. 

Even though the changing Ricci tensor strongly suggests the negative 
answer, the question is more complicated. Indeed, the intertwining oper- 
ator, constructed on ball x torus- type manifolds such that the globally de- 
fined operator is the direct sum of operators constructed on the invariant 
subspaces W 7 by the single endomorphism J 7 separately, provides an opera- 
tor intertwining the Laplacians also for the ball x ball- and sphere x ball- type 
manifolds which are in the ball x torus- type manifold. Actually, this is a 
discrete version of those constructions where one is using changing basis. 

There are many problems arising when one tries to establish the inter- 
twining of the boundary conditions for this operator on the Gordon- Wilson 
examples. First of all, one can not pass to a fixed basis and involve one- 
and two-pole functions because no fixed basis exists which is transformed 
to a well defined fixed basis by all those point transformations which define 
the intertwining operator for the subspaces W 7 . The metrics in the Gordon- 
Wilson examples are out of the touch also of the perturbation technique. 
Even though the independence theorems were guaranteed in an other way, 
they would work together just with the a deformations for which the de- 
composition v = v( a ) © holds. (This argument shows that the complete 
isospectrality can not be directly established by this discrete version of the 
intertwining operators even in case of cr-deformations, because, it is neces- 
sary to involve 1-and 2-pole functions defined by a constant basis.) In short, 
there is no way to prove the intertwining of the boundary conditions by our 
technique. 
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On the other hand, by the above considerations, a non-changing Ricci 
tensor is always a necessary condition for intertwining even the Dirichlet con- 
ditions on ball x ball- type manifolds. Therefore, these particular continuous 
operators change the Dirichlet conditions along with the Dirichlet spectra. 
Thus, they can not induce operators transforming functions defined on the 
boundaries either. In other words, operators defined by restrictions onto the 
boundary manifolds are not well defined regarding the GW-deformations. 
Even though they were introduced by suitable reductions, they would change 
the spectra also on the boundary manifolds. This phenomena strongly sug- 
gests that no other suitable operators exist and the spectra of the ballxball- 
and sphere x ball- type manifolds change during these continuous deforma- 
tions. 
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